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SLEWING ABOUT NON-ORTHOGONA L AXES 


Paul B. Davenport 
Goddard Space Flight Center 
Greenbelt, Maryland 

ABSTRACT 

Given a desired three-axis re-orientation, it is well known that there 
are (in general) twenty -four possible three-legged slews about orthogonal 
axes which will perform the re-orientation. The present paper determines 
whether the maneuver may be accomplished with three slews about arbi- 
trary (not necessarily orthogonal) axes. All of the possible three-legged 
slews are exhibited by a closed analytical procedure which requires no 
assumptions concerning the degree of non-orthogonality nor any approxi- 
mations. The classical orthogonal case (Euler angle*:) is then an immedi- 
ate consequence as a special case of the generalized solution. The de- 
velopment also includes the solutions of single-axis re-orientations by 
two successive rotations about arbitrary fixed lines. 
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SLEWING ABOUT NON-ORTHOGONAL AXES 


INTRODUCTION 

The representation of rotations via successive simple rotations about co- 
ordinate axes is probably the oldest and most frequent method employed. The 
reasors for this are threefold: (1) They are easily visualized, (2) The general 
case is obtained by three simple steps rather than a single complex one, and 
(3) They have been adequate descriptors for most problems in numerous fields. 
An arbitrary rotation can be interpreted as three simple rotations regardless 
of the actual rotation and the control or causation of rotation can be performed 
by mounting the sensors and/or actuators mutually perpendicular. Tfro-axis 
systems include office chairs, the crane, and two gimbal systems for pointing 
guns, telescopes, etc. Three-axis applications range from amusement rides to 
the precise control of ships, airplanes, and spacecraft. 

Recent requirements for precision and reliability, howe\er, have strained 
the computational luxury of traditional orthogonal axes. The imperfections of 
mounting instruments, stresses, thermal bending, etc. may cause misalignments 
of the axes. In addition to the misalignment problem, some systems are pur- 
posely designed to be non-orthogonal for reliability and/or precision, e.g., 
instruments mounted on four non-parallel faces of a regular octahedron (the 
orthogonal case corresponds to mountings on a regular hexahedron - cube). 

Here, we will take a general approach and show under what conditions 
rotations can be described by one, two, and three successive rotations about 
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fixed directed lines of arbitrary orientation. In addition, procedures are pre- 
sented for the determination of the angle of each individual rotation required to 
effect the total rotation. These angles will be referred to as generalized Euler 
angles. The development will follow what is believed to be the historical develop- 
ment of the classical Euler angles, i.e., the properties of a single axis rotation 
will be used to derive the characteristics of two successive rotations about fixed 
lines. This, in turn, forms a basis for the three axes development. In each in- 
stance the special situation of coordinate-axis rotations will be examined so that 
the results of the general case may be correlated to a familiar ov»e. 

The symbolism will be that of matrix algebra with a mixture of vector 
analysis notation. Capital English letters will be used to denote matrices. A 
vector, as used here, is a 3xl matrix and no distinction will be made in notation 
between vectors and matrices except that the latter part of the alphabet (starting 
with T) will be reserved for vectors. It is assumed that all matrices (vectors) 
are relative to some underlying orthonormal coordinate system unless otherwise 
stated. The cross product T = U x v is used in the sense that the components of 
T are formed from the components of U and V in the classical manner. Note 
that the coordinate system is not necessarily right-handed. A superscript T will 
denote transpose and the inner product of two vectors U and V expressed as 
U T V « V T U (which is equal to the dot product U • V when the coordinate system 
is orthonormal). Inverses will be described by the superscript -1. For a rota- 
tion matrix R, R" 1 = R T . 


2 



ROTATIONS ABOUV A FIXED LINE 

Given a fixed directed line defined by the unit vector U, then a rotation about 
this line is a function only of the angle, # , of rotation. The matrix of the rotation 
may be expressed as 

R ((*) - co s P I + ( 1 - co s (?) U U T - s i n 0 U, (1) 


where U denotes the skew-symmetric matrix formed from the components of 
U such that for any vector V, UV ■ U xV. I is the 3x3 identity matrix. The 
physical interpretation of eq. (1) is that a right-handed coordinate system is 
rotated (alias) by the right hand rule (when the right hand thumb points along 
the directed line the fingers point toward the positive motion) or that the vector 
is rotated (alibi) counterclockwise as viewed from the origin. The matrices of 
the rotations about the coordinate axes are given by: 



/ 


\ 


/ 





l 

0 

0 1 


' cos 

Q 

0 

- sin 6 

R, (0) = 

0 

cos # 

sin 0 

, R 2 (#) 

0 


l 

0 


1° 

-sin# 

cos 0 


^ sin 

e 

0 

cos# j 


R 3 ( 0 ) - 


cos 9 sin 0 0 

- s i n 9 cost? 0 
0 0 1 

\ 


If Y is any vector then its image Z * RY under the rotation is obtained from 
eq. (1) as 


Z = cos 9 Y + (1 - cos 9) U*YU - sin 9 V XY. 


( 2 ) 
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The following relationships are then an immediate consequence: Z • Z * Y • Y, 
Z • U ■ Y • U, and R( )U ■ R“*( ")U ■ U. The first expression is true for any 
rotation, the second holds for any rotation about U, and the last is true only for 
vectors collinear to U unless R = I (Euler's theorem). If Y and Z are given non 
zero vectors and satisfy the above conditions then the angle of rotation which 
takes Y into Z may be determined by taking the dot product of eq. (2) by Y and 
U X Y (provided U xY j* 0) respectively which gives 


co s 0 


Y • Z - (U-Y) 2 
Y- Y - ('U-Y') 2 


sin ^ 


U- (Z x Y) 

Y- Y - (U-Y) 2 


If UxY * 0 then n is arbitrary, conversely if eq. (2) holds for arbitrary 9 then 

UXY* 0. 


In order to define angles uniquely, we will make extensive use of the func- 
tion K * Tan" *(a,b) which signifies the inverse tangent of a/b where the angle 
is selected so that: (1) sgn(a) (sign of a) equals sgn(sin£ ), (2) sgn (b) * sgn(cos^), 
and (3) - v<0<>v % If a * b * 0 then \ is undefined (arbitrary). We refer to the 
angle in the above definition as the "proper quadrant" inverse tangent. Thus, 
the "smallest" rotation satisfying eq. (2) is 

6 = Tan’ 1 [U • (Z x Y), Y-Z - (U-Y) 2 ], 

all other solutions ( &) are given by e % ■ 6 + 2k v where k is an integer. 

ROTATIONS ABOUT TWO FIXED LINES 

In the preceding section it was shown that a non-zero vector Y can be ro- 
tated into a given vector Z by a rotation about a single fixed axis U if and only 
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if Z • U » Y • U and Z • Z ■ Y* Y. We now consider under what conditions a 
given non-zero vector can be rotated into another given vector via two successive 
single-axis rotations. 

Let U ( and U } be any two unit vectors such that U, x'Uj ^ 0, and 
R. ( i) -cos a I + ( 1 • cos a)lT Uj - sin a U., 

(3) 

(fi) = cos ;i I + (1 - cos . ) U j Uj - sin /3U.. 

The question to be answered then is: Given two vectors Y and Z, under what 
conditions does Z ■ R. (/^ ) R. 4 (a)Y? Since rotations do not alter length we know' 
a priori that one condition is Z • Z * Y • Y. Hereafter, we will always assume 
that this is the case. 

Now Z * (G)R.(a)Y if and only if R’ 1 ( )Z * R. (a) Y, which when expanded 

using eq. (3) yields 

( Z - U • Z U ) cos + U x Z sin .3 r IJ • Z U = 
j ) j j ) 

(Y - U. *Y U.) cos a - IJ . x Y sin a + U . • Y U , . 

v ii i ii 

By taking the dot product of both sides by U ( , IT , and U. x U. respectively we 
obtain the following three scalar equations: 


[U, - Z - (U, -11^(11. • Z}] cos /3 + U, - (U, X Z) sin P + (U. •UpfU • Z) = U. • Y, 

IT Z [U j • Y - (U. -UjXU. • Y)J cos a-U. • (U. x Y) sina + fUj •UpfU, -Y), 

(UjxU.) • Zcos P + (U. xUjWUj x 2) sin /3= (U, x U,) • Y cos a - (U, x U ) • (U, xY) sin a. 
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If we let a ■ U, • U , and 

w ( = U t *Y - a U • Z. x ( :U ( -Z-aU *Z, 

w - U • Y - a U -Y. x - U • Z - a U • Y , 

j j » i j » 

w k - (U. xU) Y. x k = (IT x U.)-Z, 

then the above equations simplify to 

x . co x is + x k sin fi - 

x. = w cos a + w, sin a 

i J ^ 

x. cos /j - x, sin = w, cos a - w sin a 

k i k j 

Note that if a * 0 the w’s and x's are merely the components of Y and Z respec- 
tively relative to the basis U. , U. , and . 

Introducing the variables «? and \ defined by 

g = w. cosa-w sin a, 

X J 

\ = Tan ’ 1 (x. , X k ), 

then equations (4a) and (4c) can be written as 

J x? + xj| sin (/^ ; X) = w. , 

Jxf + x£ COS (fi + \) = g. 

These last two equations have a common solution if xf + xjj * g 2 + w *. Thus, if 
x? + x k - wj 2 0 there exist g and /3 which satisfy both of the above equations. 
The solutions are discrete except when x 4 * x k - 0 in which case is arbitrary. 


(4a) 

(4b) 

(4c) 
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Similarly, equation** (4b) and (4c) can be expressed as 

fwj 4 w* sin ( l ♦ /*) = x , 

♦ w 2 cos (a = K> 

where ^ ■ Tan" 1 ^, .'H, )• These last two equations have a common solution if 
w ?* w k*xJ* g 2 . 

The two sets of solutions above will exist and be consistent provided that 

w 2 + w, 2 - x 2 « x? ♦ x, 2 - w 2 > 0. In this case, the equations (4) have a common 

j k j i k i — 

solution. In fact, there are at least two solutions given by: 



In terms of the proper quadrant inverse tangent 

a r Tarr 1 [(w k x j - g w^), (w j x j + g w k )] , 

£ = Tan -1 [(w. x k - g x.), (w. + g x k )]. 

The consistency of the solution above depended upon the requirement that 
w 2 + w 2 - x 2 * x 2 + x 2 - w 2 which is equivalent to the condition W T W « X T X 

j k j i k i 

where W * (w^ ,w k ) T and X * (x, ,x ; ,x k ) T . We shall now show that, indeed, 
this is the case. Obviously the condition holds when the two axes U, and Uj are 
orthonormal for then the expressions W T W and X T X represent the square of tne 
length of Y and Z respectively (assumed equal). In the general case let T and V 
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lie vectors of components of Y and Z relative to the basis U, , U, , and U, x U J# e.g., 
t, ■ U • Y, t 2 ■ • Y, and t 3 ■ (U. x U ) • Y. The above consistency condition 

is then equivalent to T T T - 2atj t ; " V^V - 2aVj v ? , where a ■ Ih • as before. 

It is possible to show that T T T - 2atj t 2 ■ Y • Y(1 - a 2 ) by constructing an ortho- 
normal basis (Gram-Schmidt Process) from U i , U , and U xU j and expressing 
the square of the length as a function of these components. A similar condition 
can be proved analogously for Z and V. This demonstrates that the condition for 
a solution to (4a) and (4c) is the same as that for (4b) and (4c) and when the solu- 
tions exist they are consistent. The above construction may also be used to 
verify that if one angle is arbitrary the other cannot be. For both angles to be 
arbitrary requires W and X to be null vectors which in turn implies that Y and 
Z be null (the components relative to the constructed basis are ail zeroK 

The condition for the existence of the solution (w|+ wj-x 2 . 0; iepends 
on the axes of rotation (U and tJ ) and verv little more can be said about existence 
until the axes (or conditions on them) are specified. If the axes are a subset of a 
complete orthonormal set (e.g., the coordinate axes) then U. • U. * 0 and the w's 
and x’s are merely permutations of the coordinates of Y and Z respectively. In 
this case, the subscripts i, j, and k can be considered as members of the set of 
integers {1 ,2,3 } and the requirement that U. and U. be distinct can be stated as 
i 4 j. There are six permutations of the three axes taken two at a time If all 
six possible choices are to be examined it is convenient to use the notation k * i'j* 
meaning k not equal to i or j. This definition plus the ’’permutation” or ’’epsilon” 
symbol e . . ( e. . ■ 0 if i ■ j, e i . « 1 if the ordered pair i, j is a cyclic per- 
mutation, and e . . * -1 if the ordered pair i,j is a non-cyclic permutation) provide 
useful tools when the rotational sequence is arbitrary. 
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l,et the ordered vectors U. , U,, U form an orthonormal triad (U ■ e . U x U 

12 3 k i j i j 

where k * i * j then all possible "two-legged" slews taking Y into Z can be expressed 
as: For i, j ■ 1,2,3 such that i 4 J let k ■ i'j', 

H = i y 2 ♦ y 2 - t 2 - t j z 2 ♦ z 2 - y 2 , 

rs ^ j J k J i k J i 

° = Tan ' 1 K'-ij Vw z , - « y,>- <y, *, * £ i, k y k >]’ 

& Tan *’ k'ti y, z k - r z i>- <y, z i 4 k z k )i. 

provided, of course, that the square root is real. 

It can be shown for the orthogonal case that at least four of the permutations 
have solutions for any Y and Z (except for the trivial case Z ■ Y) where Y • Y « Z • Z 
and depending on Y and Z there may be four, five, or six permutation solutions 
(excluding multiples of 2n there are two numerical solutions for each sequence). 

This stems from the useful fact that if the sequence (i,j) does not have a solution 
then the three sequences (j,i), ( j , k) , and (k,i) do. The non-existence implies 
z f “ y 2 > y k 2 ;> 0 and z 2 - y 2 > y 2 n 0 which in turn imply the solutions of the 
three sequences listed. Examples: (1) Y * (0,0, 1) T and Z arbitrary has only 

four permutation solutions, if z j ■ 0 then there are five, (2) Y * (1,1, 0) T and 

Z - (0,l,l) T has six (twelve numerical) solutions. 

THREE ROTATIONS ABOUT FIXED LINES 

We have now established that rotations of the form R * R^ (/SJR^a) (where 
Ri (a) and R j(/3) are rotations about fixed axes U i and U J ) are specified to within 
a positive or negative square root by the requirement that Z * RY. If one desires 
to rotate two vectors simultaneously into two other vectors by rotations about 
fLxed axes then it is clear that at least three such rotations will be required in 
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general. Since the time of Euler it has been known that relative to coordinate 
axes three such rotations (whose angles bear his name) completely define any 
rotation. In the literature, a particular (but not unique) sequence is commonly 
selected; this has caused a great deal of confusion since there are actually 
twelve possible sequences. To make matters worse some authors use left-handed 
systems while most use right-handed. Some use clockwise and others counter- 
clockwise to define the positive direction and some fail to define from where the 
motion is viewed (from the origin or looking toward the origin). Although a par- 
ticular sequence is easier to visualize (such as the yaw, pitch, and roll of a ship 
or airplane) the algebra of the general case is only slightly more complicated. 
Because of the increasing applications of non-orthogonal slewing axes and/or the 
requirement to examine all sequences (so as to select the ’’best”) we will take 
a general algebraic approach and let the reader "visualize" the results to suit 
his own application. 

It was implied above, but not explicitly stated, that two vectors and their 
images completely define a rotation. Although this is probably a widely accepted 
fact (at least in the orthogonal case) we will pursue it further for the sake of 
completeness, and since our development depends upon the result. Because of 
the orthogonality condition and the requirement that the determinant be plus one, 
the matrix of an arbitrary rotation may be written as R* (Uj ,U 2 ,U 3 ) where the 
Ui are unit vectors juxtaposed and U 3 * UjXUj. Hence, the vectors Uj and U 2 
define the rotation uniquely, they are the images of (1,0, 0) T and (0,1, 0) T respec- 
tively. Given any two vectors Y t and Y 2 such that Y j xY 2 ¥ 0 and given their 
images Zj and Z 2 under a rotation (Z. • Z ■ Y. • Y. , i,j ■ 1,2 since rotations 
preserve distance and angles) then the matrix of that rotation may be determined 
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as follows: Construct two orthonormal systems, one from the vectors Y^ and 
Y 2 , and the other from the vectors Z, and Z 2 (both of the same handedness- right 
or left). The rotation taking the first system into the second will then take Y, 
into Z i , i = 1,2. Explicitly, normalize Yj and denote this as Uj . Let W = Y 2 - 
(Yj •Y 2 /Y 1 *Y|)Y i , normalize W and denote this as U 2 . Let U 3 * U,xU 2 , and 
Rj * (U j,U 2 ,U 3 ). Derive the rotation matrix R 2 as above replacing Y with Z. . 
The rotation matrix R *= R^" 1 then satisfies Z . * RY. , i = 1,2. The uniqueness 
of R follows from Euler's theorem (a nontrivial rotation has only one real 
eigenvector). 

From the discussion above we have ascertained that the practical question - 
Under what conditions can two vectors, a coordinate system, or a rigid body be 
rotated to some desired orientation via three fixed-axis rotations ? - can be 
formulated algebraically as follows: Under what conditions can an arbitrary ro- 
tational matrix R be represented as R ■ R n ( ^)R(> (<£) where R . (i * ( ,m,n) 
denotes a rotation about a fixed unit vector U i ? It is assumed that m is different 
from n and f (otherwise the product degenerates into only two fixed-axis rotations) 
but that n and P may be equal. 

Assume R = R n (0)R m (^)Rp(0) with the conditions above, then since R UJU, 

= Rj 1 (V)U, = Uj (i <* f ,m,n and \ = 4 > , we obtain the following results: 


R U e = R n (./>) R m (?) U f , (5a) 

R' 1 U n = Rf' (0) R‘> (?) U n , (5b) 

R Rj 1 U) U. = R n (<p) U m , (5c) 

Rn R „ W U m = Rj ! (4>) U m . < 5d > 
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Each of these equations has the property of eliminating one of the three gener- 
alized Euler angles. If the U ( are coordinate axes then the left side of eq. (5a) 
is merely the P th (p » 1,2,3) column of R and the left side of eq. (5b) is the n th 
row of R represented as a column vector (the n th column of R T ). 

Conversely, assume that equations (5a) and (5b) hold and set Yj * Up , Z 2 ■ U n , 
Z 1 *= R (^RJ )U f , and Y 2 = Rj 1 <*)R ' nl ‘(f , )U n . Then Zj-Zj-Yj.Yj-Zj.Z,- 
Y, • Y, = 1 and Z, • Z 2 = Y, • Y 2 - U T R m (0)U f . Therefore, If \ xY 2 4 0 
(R m ( >? )Uf> 4 ±U n ) it follows (by a previous argument) that the R satisfying Z. * RY. 
(i = 1,2) (equivalent to equations (5a) and (5b)) is unique. This means that the 
simultaneous solutions of equations (5a) and (5b) provide the solution to R - R n 
(v )R m (‘ ; )Rp(^) except in the isolated cases where R m (^) Uf is collinear with U n 
(e.g., when ( - n and 0 * 0). This is indeed felicitous since both equations are of 
the same general form which was solved in the previous section. We merely 
have to demonstrate that the two solutions are consistent under the present 
assumptions. 

Given Uj , U 2 , U 3 , R, and a particular rotational sequence, then the solution 
of eq. (5a) is obtained directly by the techniques of the previous section by letting 
Y * Up , Z * RUp , i = m, j = n, 6 » a, and 0 . Eq. (5b) is solved similarly. 
Before doing this, however, it is convenient to introduce some intermediate 
quantities which occur often and which help to simplify the final result. 

Let 

V 1 = U 2 x U 3- V 2 = U 3 * U r V 3 = U , 

flij = a.. = U t -Uj, (i, j = 1, 2, 3) 

f=U • V =U • V = u • V 
U 1 V 1 U 2 V 2 u 3 v 3’ 
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and when R is obtained 


•>, , = U, • (R Uj ) = U 1 -(R-' Up. 
c l( = V ( -(R U ^ (i. j = 1. 2, 3) 

d. . = U, • (R Vp = \ l - (R _1 Up 

Define the matrices M * (U ,U ,U ) and N ■ (V ,V ,V ) (juxtaposition of vectors) 

i • *3 i A J 

then the above definitions may be expressed as: 

A = M T M, B = M T R M, C = N T R M, I) = M T R V. 

The matrices C and D may also be written in terms of A and B, i.e., 

C = d (M) A -1 B, D = d (M) BA' 1 , 

where d(M) denotes the determinant of M. If the Uj are the coordinate axes then 
A=N=M=I, B=C*D*R, and f * 1. 


Returning to the solution of eq. (5a), the proper substitutions give 


w 

= a n 

- a 

b n t 

X 

= b n 

- a 

b n 

m 

m t 

mn 

nr 

m 

mr 

mn 

nt 

W 

= a Q 

- a 

a o , 

X 

= b rt 

- a 

a 0 

n 

ni 

mn 

r m 

n 

nt 

mn 

i m 

W k 

= e nf 

f, 


X k 

= £ 

mn 

c kf’ 



where k = m’n'. If 

h = w 2 + w 2 - x 2 = x 2 + x 2 - w 2 > 0, 

n k n m k m ’ 

(the equality always holds) then g = ± J”h and 

6 = Tan' 1 [(w k x n - g w n ), (w n x n + g w k )] , 

4 , = Tan’ 1 [(w m x k - g xj, (w m x m + g x k )]. 


( 6 ) 


( 7 ) 


(8a) 

(8b) 
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Also, 


g - w. cos ^ - w sin 

" k n 


x, cos - x s i n 0. 


, let 

Y = 

u 

, Z 

- r m u , 

i ■ m, 

j- 

f, 

a 3 



n 


n 




w. = 

i 

3 mn 

- 


b „f 

x - 

i 

1) 

nm 

n p 

w f = 


- 


a - w , 

mn n 

Xf * 

b 

n 

f - 

a ( 

w = 

€ fn 

f 

= - 

V 

X - 
p 

C 

m 

f d 

np’ 


(9) 


where p » Pm’ (the subscript i is for identification only having no numerical sig- 
nificance, i.e., w, 4 w m ). The condition for a solution is 


w- 


w 2 - x 2 = x? 4. x 2 - w? z 0 


(10) 


f T " _ "f ‘ p 

(again the equality is assured by the initial conditions), but the first expression 

is identical to h defined above in the solution of (5a). Thus, if (5a) has a solution 
so does (5b) and vice versa. Furthermore, the values of 0 from (5b) can be 
brought into agreement with those already obtained from (5a) by proper selec- 
tion of the sign of the square root of h. This gives 


g' = w cos d - w p sin (- 8), 

= x p cos 0 - x, sin ( - 0) , 

- - (w k cos 9 - w p sin 9) = - g, 
and = x . cos 0 + x sin (- 0) - w . 

I P \ t - * i 

It then follows that 


4 > = Tan" 1 [- (w. x p + g x.), (w. x. - g x p )J . 


( 11 ) 
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In summary, if h, as defined above, is non-negative then (5a) and (5b) have a 
common solution given by the above equations. This implies that R can be fac- 
tored as R = R„(, )R m ( )Rj>(sO provided R yn ( , )U ( > ^ ±U n . 

At least one of the arguments of the proper quadrant inverse tangent in (8a) 
is different from zero (either w or w k not zero) which means the generalized 
Euler angle always has well-defined discrete values when the solution exists. 
If this solution is such that R t ( )Up = ±U n then (5a) and (5b) both reduce to 
RU } > = ±U n regardless of > or </•. Conversely, if RUp * ±U n then (5a) and (5b) 
degenerate into the single condition R m (^)Up * ±U n . Hence, we have ascertained 
that the common solution of (5a) and (5b) gives the desired factorizations of R 
except when RUj> = ±U n in which case \p and / are arbitrary. In general, there 
are two discrete solutions (neglecting multiples of 2v) corresponding to the two 
square roots of h. 


When the axes of rotation are the coordinate axes then the expressions above 
can be greatly simplified. In this special case, 


h - 


P 


r kT 


r 2 4. r 2 

ntn np 


= 1 - 


'n 2 P 


t nI h , 


* = Tan ‘ 1 r „p>’ 

" Ta,r ‘ 'K ( r nf - *nf . ( i> n f + <„( g) . 
'i ' r Tan ‘' r m f 4 f .n r ke^ 


k = m f n\ p = f , m* , and b 0 is the Kronecker Delta ( * 1 if n * P and & 0 * 0 

r n l ni nt 

if n 4 P). There are twelve possible sequences (six for n = P and six where 
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n ■ Pm') and each has two solutions except when h ■ 0. In this latter case, / 
and 0 are undefined and "■0ifn«for ,, *± ”/2 if n ■ Pm'. This corresponds 
to the condition RU|> * ±U n . 

The degenerate case (RUj, * ±U n ) can be handled by the use of (5c). In fact, 
it can be demonstrated that the simultaneous solution of (5b) and (5c) provides 
a solution to the factorization of R into three fixed-axis rotations whenever a 
solution exists. The argument is analogous to that used for (5a) and (5b); the 
consequence that R - R r (0)R m ( ' )Rp(^) (5b) and (5c) hold requires only that 

U n 4 ±U m which was assumed to be the case. Eq. (5c) provides the angle i> as 
a function of <fi which can be evaluated after <p is determined from (5b) or 
assigned an arbitrary value when (5b) is degenerate. This functional relation- 
ship may be derived by employing familiar techniques; expand (5c) using the 
expression for a rotation about a line and then take appropriate dot products of 
each side. The result is: 

O - a nm'> sin * = - < V n * U «> ‘ K Z m " 8 f m Zf) COS 0 + Zf * Z,, sill 0 + Zf] , 

(1 - O cos^ = U m - l(Z m .„, m Z p ) cos 0 + Z f * Z, sin 0 + a Pm Zf - a nm UJ , 

where Z x = RU A and a, . = IT • U j as before. Introducing the matrix E * N T RN 
(e. . = V. • RV. ) where the V. and N were defined previously, the relationship 
between $> and 0 may be described as follows: 

0 = Tan' 1 [y (0), z (0)], (12) 

where 

y W = £ .n ^ C km - a f* COS <t> * e f. e kp sin * + a e™ c k6 ] - 

Z - a u b «£> COS * * d mp sin * ♦ a f„ b .f * a ’ m - 
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rind all other quantities have been defined heretofore. The transpose of the 
matrix E is equal to the adjoint of the B matrix. In the elassical situation of 
rotations alxmt coordinate axes this last formula reduces to 


Tn:i 


( r, cos 

tr ii km 


*'r 


r. sin ;), (r cos 

k p tn m 


« c r 

• m nip 


s i n 




Although the derivation, including existence proofs, of the formulas con- 
tained herein has been somewhat lengthy and tedious, the application of the re- 
sults is straightforward. The matrices M,N, and A are constant depending only 
on the fixed axes of rotation (if only a single fixed sequence is being considered 
many elements need not be evaluated). The matrices 3, C, D, and E are easily 
determined by matrix multiplication once R is known. The inequality (7) or (10) 
establishes the existence or non-existence of the solution once the parameters 
in (6) have been determined. Eq. (8a) provides the well defined angles for the 
middle rotation. The angles for the first rotation are evaluated from (11) then 
the angles for the last rotation may be obtained by (12) or alternately by (8b) if 
the first angle is well defined. 
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